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The periodized cluster dynamical mean-field theory (PCDMFT) combined with exact diagonaliza- 
tion as impurity solver has been applied to the half-filled standard Hubbard model on the honeycomb 
lattice. A correlated Dirac semimetal is found for weak interactions and it transforms into an an¬ 
tiferromagnetic insulating phase for strong interactions via a first-order quantum phase transition, 
not intervened by a spin liquid phase in between. In this application, the PCDMFT introduces the 
partial translation symmetry, but cures well the problem due to the translation symmetry breaking 
in the cluster dynamical mean-field theory studies for the same model, which give rise to a spurious 
insulating phase in the weakly interacting region. 
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I. INTRODUCTION 

In past decades, much effort has been devoted to 
searching for materials and models that show a quan¬ 
tum spin liquid, in which the local moments fluctuate 
even at zero temperature. Such a state is a genuine 
Mott insulator without spontaneously breaking spatial 
and spin symmetries and can show remarkable emergent 
phenomena such as non-trivial topological order, frac¬ 
tionalized charges and excitations. A common idea is 
that geometric frustration may prevent orders and lead 
to spin liquids 1 , such as in models on kagome lattices. 
Another proposal is that a spin liquid may emerge nearby 
a Mott transition 2,3 , such as in the Hubbard models on 
triangular or honeycomb lattices. Meng and coworkers 
4 attempted the later idea. They performed large-scale 
quantum Monte Carlo (QMC) calculations on a standard 
Hubbard model at half-filling on the honeycomb lattice 
at zero temperature. After the obtained QMC results for 
finite-size clusters containing up to 18 x 18 x 2 sites were 
carefully extrapolated to the thermodynamic limit, they 
found a spin liquid state between the Dirac semimetal at 
weak interactions and the antiferromagnetic insulator at 
strong interactions. But, later this result was challenged 
by Sorella et al. 5 with the similar QMC calculations on 
larger clusters with up to 36 x 36 x 2 sites and by Assaad 
and Herbut 28 with an advanced QMC technique. 

On the other hand, the existence of spin liquid on 
the honeycomb lattice has also been studied 6 ” 11 by using 
quantum cluster methods 12 , which map a lattice model, 
a Hubbard model for example, onto a cluster of impu¬ 
rity sites coupled to a set of bath sites determined self- 
consistent ly or by a variational principle. Wu et al. 6 
used cluster dynamical mean-field theory 13 (CDMFT) 
combined with continuous-time quantum Monte Carlo 14 
(CTQMC) as impurity solver, while Liebsch 7 also used 
CDMFT but with exact diagonalization 15 (ED) as impu¬ 
rity solver. They both did calculations at finite tempera¬ 
tures and showed the existence of a spin liquid and a Mott 
gap at intermediate interactions. However, performing 
also ED-CDMFT but at zero temperature and with elab¬ 


orate numerical analytic continuation on the Matsubara 
Green’s functions, He and Lu° found that the existence 
of the single-particle spectral gap extends actually to the 
weak interaction limit, i.e., for all interactions U > 0. 
Later, this result was reproduced by Seki and Ohta 10 , 
who used variational cluster approximation 16,17 (VCA) 
with also ED as impurity solver. Since it is known that 
the Dirac semimetal is stable for weak interactions 18,19 , 
the validity of the application of CDMFT-like methods 
to the present model is questioned 11,20 ” 22 . 

In all works using quantum cluster methods mentioned 
above, the impurity clusters consist of six sites forming 
a ring with the same rotation symmetries as the hon¬ 
eycomb lattice. When the ED is used as the impurity 
solver, only six bath sites couple to the impurity cluster, 
which is close to the ED’s computational limit. Apply¬ 
ing VCA, CDMFT, and the cluster dynamical impurity 
approximation 23 (CDIA) at zero temperature with the 
ED as the impurity solver to honeycomb and square lat¬ 
tices, Hassan and Senechal 11 argued that one bath site 
per boundary impurity site in ED CDMFT and ED VCA 
is inadequate and leads to the absence of the semimetallic 
phase for the honeycomb lattice. In contrast, when they 
adopted two- or four-site impurity clusters with two bath 
sites per boundary impurity site, a first-order transition 
at a finite interaction from the semimetal to the antifer¬ 
romagnetic insulator was obtained. However, this argu¬ 
ment was refuted by Liebsch et al. 20,21 , who showed that 
CTQMC CDMFT, which uses a continuum of bath sites, 
yields the same results as the ED CDMFT at finite but 
low temperatures for the honeycomb lattice. Moreover, 
with the detailed analysis 20 ’ 21 on the ED CDMFT cal¬ 
culations, they attributed the absence of the semimetal¬ 
lic phase to the translation symmetry breaking in the 
CDMFT. 

In order to avoid the translation symmetry breaking, 
Liebsch and Wu 20 resorted to dynamic cluster approxi¬ 
mation 24,25 (DCA), which is a quantum cluster method 
preserving translation symmetry naturally. After using 
the DCA, the semimetallic phase at weak interactions 
was successfully recovered. However, they found that 
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this semimetallic phase is still robust near U = 6t at 
which the CDMFT and large-scale QMC calculations al¬ 
ready give rise to a Mott insulating phase. This shows 
that the DCA overemphasizes the semimetallic behav¬ 
ior, which is mostly because its condition on ensuring 
translation symmetry is too rigid for the description of 
correlations within the impurity cluster. Thus neither 
DCA nor CDMFT can describe the correlations properly 
in the whole interaction range on the honeycomb lattice. 

Here we revisit the problem, namely, the quan¬ 
tum phases of the half-filled standard Hubbard model 
on the honeycomb lattice, with periodized CDMFT 
(PCDMFT) 26 ’ 27 , to realize a proper description of the 
semimetallic phase at weak interactions and meanwhile 
retain CDMFT’s advantage in well describing correla¬ 
tions at strong interactions in the framework of quan¬ 
tum cluster method. The PCDMFT is a variation of 
CDMFT partially restoring translation symmetry. As 
shown below, the PCDMFT indeed well reproduces the 
semimetallic phase at weak interactions, separated from 
the antiferromagnetic insulating phase at strong interac¬ 
tions by a first-order phase transition with the U c ~ 3.7 1, 
and there isn’t a spin liquid emerging near the transi¬ 
tion. These are in agreement with the results from the 
large-scale QMC calculations by Sorella et al. 5 

II. MODEL AND METHOD 


(a) (b) 

FIG. 1. (Color online) (a) The honeycomb lattice is a lat¬ 
tice with each primitive unitcell consisting of two sites. The 
CDMFT/PCDMFT considers it as a superlattice with six-site 
rings as supercells, for example, (b) The CDMFT/PCDMFT 
maps the Hubbard model on the honeycomb lattice onto an 
impurity model by considering a supercell as an impurity clus¬ 
ter coupling to a non-interacting bath (schematically shown 
by the shaded area) which approximates all other supercells. 

As in those works discussed above, we study the half- 
filled Hubbard model on the honeycomb lattice as follows, 

H=-t^2 (cJo-Cjcr + c] a C ia ) +UJ2 n rt n ih (!) 

* 

where a =j" or Ci G annihilates an electron with spin 
a at the ith lattice site, is the number of electrons 


with spin a at site i, and ( ij ) denotes a pair of nearest- 
neighbor sites i and j. The studied honeycomb lattice is 
shown in Fig. 1. 

The CDMFT considers the lattice as a superlattice 
composed of supercells, each of which consists of a cluster 
of sites, as shown in Fig. 1. Then the CDMFT takes a 
cluster as the reference cluster, or impurity cluster, and 
approximates all other clusters by a non-interacting bath 
coupling to the impurity cluster, i.e., maps the original 
problem onto an impurity problem. The bath is deter¬ 
mined self-consistently. The CDMFT neglects nonlocal 
self-energies between different clusters, i.e. E MI/ = 0 for 
H v where ^ and v enumerate clusters, and further 
identifies the impurity cluster self-energy E 7 in the im¬ 
purity problem with the local self-energy E M/i in the orig¬ 
inal problem, i.e., = E 7 . It is this assumption that 

makes the CDMFT break translation symmetry. 

To establish a complete correspondence between the 
lattice problem and the impurity problem, the CDMFT 
introduces a Weiss field Q defined by 

I ?" 1 = ( G ^)" 1 + E w , ( 2 ) 

where E MM and G^ are actually independent of /j, be¬ 
cause of the translation symmetry of the superlattice. 
And G w can be derived from 

G(K, iuj n ) = [iuj n - H 0 (K) - E(K, *w n )] _1 (3) 

with K the superlattice momentum, i7o(K) the hopping 
coefficients of the model on the superlattice in the K- 
space, and E(K) = E 7 independent of K by a simple 
derivation. Finally, with the Dyson equation of the im¬ 
purity model E 7 = (Gg) _1 — (G 7 ) -1 , the self-consistent 
equations of the CDMFT are closed by identifying the 
non-interacting Green’s function of the impurity model 
Gq with the Weiss field Q, i.e., 

(Go) -1 = Q 1 ■ (4) 

In numerical realization of the CDMFT, the self- 
consistent equations are solved by an iterative algorithm. 
One begins with a guessed E, for example, the one ob¬ 
tained from other calculations or simply by setting E = 0. 
Substituting it into Eq. (3) and Eq. (2), one can obtain 
G and Q, successively, and then Gq from Eq. (4). Mean¬ 
while, the impurity problem is independently solved by 
using an impurity solver such as ED or QMC to yield the 
Green’s function of the impurity sites G 7 . Thus, accord¬ 
ing to the Dyson equation of the impurity model, one 
obtains E 7 , which then gives a new E for the original 
lattice model. This process will be iterated until conver¬ 
gence. 

When the ED is used as an impurity solver, the bath 
is discretized, i.e., described by a set of non-interacting 
sites coupling to the impurity sites. In a general form, 
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the impurity model reads 

11 ' = E tijcLcjo- + U E n if n 4 

i,jG cluster, cr iGcluster 

+ ^2 £ kc{c k + ^2 (Vifccjcfc + H.c), 

fcGbath iGcluster,fcGbath 

(5) 

where and 17 coincide with those in the original lattice 
model. And £k and Vjfc enter into Gg with 

(Gg) _1 (iWn) = *w n - t c - V(iu n - s)~ 1 V\ (6) 

where e is a diagonal matrix with diagonal elements 
£fc, while the elements of matrix V and t c are V)fc 
and tfj, respectively. To determine £k and Vik, we fit 
(G'g)- 1 (iw n ) with Q~ l (iujn): he., minimize the distance 
between (Gq) - 1 and t/” 1 by adjusting unknown param¬ 
eters £fc and Vifc. 9,15 

Now we turn to the PCDMFT. As introduced above, 
the CDMFT breaks translation symmetry. By inspect¬ 
ing the self-consistent equations, we see that the lattice 
quantities £(K) and G(K) violate the translation sym¬ 
metry of the original lattice because of directly setting 
equation £(K) = £ 7 . To restore the translation symme¬ 
try, in the PCDMFT we construct the lattice self-energy 
£ through periodizing £ 7 as follows, 26,27 

E ( k ) = ^E e “ lkR “ E V kR3 > (7) 

aft 

where a (/3) enumerates primitive unitcells in the cluster, 
R q (Rg) is the lattice vector for primitive unitcell a (/3), 
k is the lattice momentum, and N c denotes the number 
of primitive unitcells in the cluster. Accordingly, Eq. (3) 
changes into 

G(k, ito n ) = [iu) n -H 0 ( k) - £(k, iw n )] _1 . (8) 

And in comparison with the Dyson equation of the im¬ 
purity model, Eq. (2) is changed into 

G- 1 = (G^)- 1 + £ 7 , (9) 

because £ w = £ 7 no longer holds after periodizing the 
self-energy. From now Eq. (9) will work with the Dyson 
equation of the impurity model to guarantee the iterative 
consistency. 

Following the previous studies 6-10,20 in the literature, 
we choose a six-site ring as the impurity cluster for the 
honeycomb lattice, as shown in Fig. 1. However, this 
choice is troublesome for the periodization (Eq. (7)) since 
such an impurity cluster cannot be divided into a set of 
primitive unitcells. To overcome this difficulty, we first 
decompose Eq. (7) formally into 

4 = A £ Ke, (10) 

C -R Q =R/y 

£(k) = ^£^e ik - R G (ll) 


and further generalize them into 

= W E ( 12 ) 

C R/3—R Q =R-y 

£ afe (k )=J2Kab^^, (13) 

7 

where a and b enumerate sites in primitive unitcells; es¬ 
pecially, in a and b denote the sites in unitcells a 

and /?, respectively. When the impurity cluster can be 
divided into a set of primitive unitcells, Eqs. (12) and 
(13) are equivalent to Eqs. (10) and (11), respectively. 
However, when the impurity cluster not, the latter no 
longer works while the former still works. In this case, 
we still divide the impurity cluster into a set of primitive 
unitcells, in which some unitcells are partly included. We 
then assign zero to the impurity self energy with a site 
beyond the impurity cluster, represented by Eq. (12). 


III. RESULTS 

We have applied the PCDMFT to the half-filled stan¬ 
dard Hubbard model on the honeycomb lattice. To study 
the ground state properties, we used the ED as an im¬ 
purity solver. Thus, we calculated the same physical 
quantities as those by the CDMFT reported in Ref. 9, 
and then compared them between the PCDMFT and the 
CDMFT. Specifically, the calculated quantities are the 
density of states, single-particle gap, staggered magneti¬ 
zation, and double occupancy, respectively shown in the 
following Figs. 2, 3, 4, and 5. 



co/t 

FIG. 2. (Color online) Density of states (DOS) calculated 
by the PCDMFT. The linear behavior near the Fermi energy 
shows a correlated Dirac semimetal. 

For small interactions U, the density of states (Fig. 2) 
shows a perfect linear behavior near the Fermi en¬ 
ergy, which means that the expected correlated Dirac 
semimetal is successfully produced. Correspondingly, 
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FIG. 3. (Color online) Single-particle gap A/t calculated by 
the PCDMFT. It shows hysteresis behavior around U c = 3.2— 
3.7 t, indicating a first-order quantum phase transition. For 
U > U c , the system opens up a gap, indicating an insulating 
phase. 

both the single-particle gap A/t (Fig. 3) and the stag¬ 
gered magnetization m (Fig. 4) are zero. In comparison, 
the CDMFT results in Ref. 9 show a nonzero A /t for any 
finite interaction. 

As the interaction U increases, the system encoun¬ 
ters a first-order phase transition from the semimetallic 
phase to an antiferromagnetic insulating phase. The two 
phases coexist in the transition region of U c i < U < U c i 
with U c \ = 3.2 1 and U C 2 = 3.7 1. And the characteris¬ 
tic hysteresis behavior for a first-order phase transition 
shows up in the staggered magnetization (Fig. 4), the 
single-particle gap A/t (Fig. 3), and the double occu¬ 
pancy D (Fig. 5) as functions of U, respectively. The 
converged result depends on the initial condition (the 
lattice self-energy E in the present study) for the self- 
consistent iterations. Specifically, the converged result 
is the semimetallic phase if the initial E is equal to the 
converged E from a previous calculation with a small 
interaction or is simply set to zero; in contrast, the con¬ 
verged result is the antiferromagnetic insulating phase if 
the initial E is equal to the converged E from a previ¬ 
ous calculation with a large interaction, as shown by the 
arrows in the figures. 

For U > U c 2 , the converged result no longer depends 
on the initial E and is always the antiferromagnetic in¬ 
sulating phase, indicated by both the nonzero staggered 
magnetization (Fig. 4) and the nonzero single-particle 
gap A/t (Fig. 3). This is consistent with the result by 
the CDMFT in Ref. 9. 

As seen from the above, the main difference be¬ 
tween the PCDMFT and CDMFT results 9 is that the 
PCDMFT produces a correlated Dirac semimetal at 
small U, while the CDMFT produces an insulator for 
any finite U. The other results are similar, such as the 
double occupancy D for all interactions. Nevertheless, 



FIG. 4. (Color online) Staggered magnetization calculated 
by the PCDMFT. Hysteresis behavior is shown around U c = 
3.2 — 3.7 1, indicating a first-order quantum phase transition. 
For U > U c , the staggered magnetization is nonzero, indicat¬ 
ing an antiferromagnetic phase. 


U c is about 3.7t for the PCDMFT, which is remarkably 
closer to the values of U c = 3.8 1 5 and U c = 3.7 t 2S ob¬ 
tained from the large-scale QMC simulations, compared 
with that U c ~ 4.6 1 for the CDMFT 9 and U c ss 3.3f for 
the DC A 20 . 



FIG. 5. (Color online) Double occupancy calculated by the 
PCDMFT 


IV. DISCUSSION AND CONCLUSION 

The cluster dynamical mean-field theory (CDMFT) 
breaks the lattice translation symmetry by construction, 
which leads to a spurious phase that an insulating phase 
emerges at any small but finite interaction for the half- 
filled Hubbard model on the honeycomb lattice. To solve 
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this issue, we have employed the periodized CDMFT 
(PCDMFT) which improves the CDMFT to partly re¬ 
store the translation symmetry by periodizing the impu¬ 
rity cluster self-energy. It turns out that the PCDMFT 
successfully clears up the spurious insulating phase at 
small interaction U and gives rise to a correlated Dirac 
semimetal instead. As U increases, the system takes 
a first-order phase transition with U c ~ 3.7 1 from the 
semimetal to an antiferromagnetic insulator. And there 
is no a quantum spin liquid emerging near the transi¬ 
tion. Basically, these results are in good agreement with 
the large-scale quantum Monte Carlo (QMC) results by 
Sorella et al. 5 

The PCDMFT inherits the CDMFT but meanwhile 
keeps translation symmetry like the dynamic cluster 
approximation (DCA). Thus it may be considered as 
a method interpolating between the CDMFT and the 
DCA. For the present model study, the DCA 20 underes¬ 
timates the transition point U c while the CDMFT 9 over¬ 
estimates it. In comparison, the PCDMFT gives the one 
most close to the large-scale QMC result 5,28 . This indi¬ 
cates that the PCDMFT makes a better balance between 


the descriptions of weak and strong correlation effects 
on the honeycomb lattice than the other two. Although 
these quantum cluster methods all guarantee to predict 
the same result by using an enough large impurity clus¬ 
ter, they may differ with each other significantly in pre¬ 
diction by using a small impurity cluster and thus should 
be used complementarily in practical studies. 
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